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1 . Introduction 


Galvano- and thermomagnetic phenomena in conductors are 
usually described with the help of kinetic coefficients which are 
determined with the following expressions for electric field E 
and the electron fraction of the heat flow q [ 1 ] : 

[ kx j ]+c*VT+QB[ kxVT ] , ( 1 ) 

q =n j+E/cB[ kx j ]-/cVT+LB[ kxVT ] . (2) 

Here p xx and p xy are com P onents of the resistivity tensor; j is a 

current density; k=B/B is a unit vector parallel to magnetic 
field B; a, Q, IT, E, k , and L are Seebeck, Nernst, Peltier, 
Ett ingshausen , heat conductivity and Righi-Leduc coefficients, 
respectively. In accordance with the Onsager principle, some of 
these coefficients are interconnected: 

n=aT, E/c = QT . (3) 

On the other hand, there is also a connection between the 
electric field and the heat flow in the superconductor mixed 
state. If one assumes that both the electric field and the heat 
flow arise only under the fluxoid motion then we have Josephson 
formula for the electric field: 

E =-|[v f xB] (4) 

and Huebener relationship for the heat flow [2]: 

3 =£ > Tn f V f - (5) 

where v f , n f =B/$> 0 are an average motion velocity and a fluxoids 

density (3> 0 =hc/2e is a magnetic flux quantum) , and is an 

entropy of a unit section of a fluxoid. It follows from Eqs.(4) 
and (5) : 


q=T(cS^/* 0 )[Exk]. 


( 6 ) 


318 



For given values of the magnetic field and temperature the 
velocity of the fluxoid motion v f depends only on the current 

density j and the temperature gradient VT. For a general case the 
dependence of v f on j and VT has the form: 

v f =a a j + a 2 [jxk] + a 3 VT + a 4 [VTxk] (7) 

(a specific form of coefficients a 3 , a 2 , a 3 , and a 4 depends on 

the choice of the fluxoid motion model, see below) . Substituting 
Eq.(7) into (4) and (5), we obtain 

E=(a 2 B/c)j + (ajB/cMkxj] + a 4 BVT - a 3 B[kxVT], (8) 

q=(n f TS^a 1 )j - (n f TS^a 2 ) [ kx j ] + a 3 VT - aJkxVT]. (9) 

The comparison of these relations with general expressions (1) 
and (2) for E and q and the use of the Onsager equation (3) lead 
us to 


“s^VV^y • % B = -< C VV4 • < 10 > 

where index s marks kinetic coefficients relating to a 
superconductor in the mixed state. Eqs.(10) which are a 
consequence of the connection (6) between q and E in the mixed 
state of a superconductor forecast a certain correlation between 

s s 

"heat" and "current" coefficients (a , Q and p , p ) and give 

S s xx x y 

a possibility of checking the correctness of the conception 
involved . 

The most detailed experimental researches of kinetic 

coefficients in the superconductor mixed state have been carried 

out for high temperature superconductors [3,4] where their value 

is essentially higher and the temperature range accessible for 

measurements is noticeably wider than those for traditional 

superconductors. The researches revealed the real connections 

between kinetic coefficients to be strongly different from those 

predicted by relations (10). Thus, for instance, Hall resistance 

p in the mixed state is often a nonmonotonic and even 
r xy 

sign- reversing function of temperature and magnetic field. 
Meanwhile Seebeck coefficient (thermal power) usually rises 
monotonically with the T or B increase and its behavior, as a 
rule, resembles that of the longitudinal resistance :q a p^. 

Besides, according to Eq.(10), the Seebeck coefficient goes to a 
zero on the normal state transition (since in this case entropy 
goes to a zero) , although this coefficient obviously should 

take the value peculiar for the normal state. All this testifies 
to the fact that Eq.(10) is based on incorrect premises. It is 
clear that these are only Eqs . (4) and (5) (ensuing from an 

assumption that the electric field generation and the heat 
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transfer and, consequently, all the kinetic phenomena in the 
superconductor mixed state are related solely to the fluxoid 
motion) may be erroneous. 

This approach is applicable in full measure only far from the 
superconducting transition, where fraction N =n /N of normal 

n n 

quasiparticles is small (N and n n are normal quasiparticles 

concentrations at T>T c and TcT, , accordingly). In conformity with 

the weak coupling BCS theory such particles originate only due to 
their excitation via energy gap A, and hence in low temperature 
superconductors (with a weak electron- phonon interaction) normal 
quasiparticles appear in a noticeable quantity only in a direct 
proximity to T c when the width of superconducting gap A becomes 

close to (or smaller than) thermal energy kT. Thus, this approach 
is valid almost everywhere but for a narrow temperature range 
near the transition. 

A different situation occurs in high- temperature super- 
conductors with a marked tendency towards a strong 
electron- phonon interaction resulting in the finite density of 
states within the energy gap (5). In this case normal 
quasiparticles arise mainly owing to filling these states and 
their concentration is determined by a power- law dependence 

(e.g., by law N n s«(T/T c ) 4 , known for a Gorter-Casimir two- liquid 

model [6]) rather than by temperature exponential dependence 
N n (T) a exp(-A/kT) . The temperature range, where fraction of 

normal quasiparticles is large, grows much broader and there is 
no groundj to neglect their contribution into different kinetic 
phenomena . In Ref . [ 5 ] numerous references are made to the works 
which point out that the temperature dependences of thermal 
properties of high temperature superconductors near T c are much 

better described by the Gorter-Casimir phenomenological model 
than by weak coupling BCS theory. In this connection the 
description of kinetic properties of the mixed state of high 
temperature superconductors requires a simultaneous and 
self-consisted account for the motion of fluxoids and normal 
quasiparticles . 

The most straightforward way here is to proceed from an 
assumption that corresponding contributions into kinetic 
coefficients are independent and additive. And it exactly this 
way which is usually taken (see, e.g., Refs. [7 and 8]). This 
approach presupposes, in particular, that the presence of normal 
quasiparticles effects in no way the fluxoid motion (and, 
consequently, the electric field in a superconductor) . This key 


Here by the term "normal quasiparticles" we understand not the 

particles which are "bound" in the fluxoid cores (with the 

concentration of N) but those which are outside the fluxoid cores 

and are "free" (their concentration is n <N) . 

n 
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assumption is corroborated with no due arguments and is simply 
postulated although it does not look absolutely obvious. So one 
needs to get rid of this limitation and to construct a 
self-consisted theory of galvano- and thermomagnetic phenomena in 
the superconductor mixed state within a frame work of a single 
model. Since at present the fluxoid motion is described 
approximately even in neglect for normal quasiparticles [9] a 
phenomenological approach analogous to that used in the classical 
work by Gorter-Casimir [6] might be used. 

2. Seebeck and Nernst Effects 


Average density j of the normal quasiparticle current is 
determined by the standard relationship: 

j =a N E+a N [ Exk]-/? . VT-/3 [ VTxk] , (11) 

J n xx xy L J r Nl r H2 

N N 

which takes into account that conductivity (a^ x and o ) of a 

system of normal quasiparticles (metal) is proportional to their 
concentration N n while Seebeck (a N ) and Nernst ( Q N ) coefficients 

2 N N 

are independent on it . Here ^Ni =Q N°xx +< ^N Bo xy an< * 

B =a ct N -Q Bo N . The electric field in a conductor with finite 
r N2 N xy xx 

conductance is determined by well-known relations, namely, the 
Ohm's law 


E = p j +p [kxj ]+q VT+Q B[ kxVT ] 

r xx J n r xy l J n N N 


( 12 ) 


for a normal metal and Josephson’s formula (4) for a 
superconductor in the mixed state without normal quasiparticles. 

Here p N and p N are components of resistivity tensor of a normal 
r xx r xy 

metal. How could one generalize Eqs.(4), (12) corresponding to 

two limiting situations ( N n = 0 and N n = l) ^ or a common case when 

0<N <1? In the absence of an exact answer to this question one 
n 

may try the interpolation formula 


E =- f s {5 ( V B ] 

where weight 


1+f (p N j +p N [kxj )+o 7T+0 B[ kx VT ] i 
J n [ r xx J n r xy l J n N N 

factors f and f have to satisfy 

s n 


(13) 


conditions 


p 

In normal metal a, Q cc kT / c [1] where Fermi energy c p is determ 
the total concentration of charge carriers and is independent on N n . 
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f (N =0) =1, f (N =0)=0, f (N =1) =1 which ensure the transfer of 

Eq . (13) into "correct" expressions (4) and (12) in the above 
cited limiting cases. The expression for an electric field which 
is obtained after substitution of Eq.(ll) into (13) contains no 
terms proportional to VT (it is required for obvious equality E=0 
to take place at v f = 0 ) and has the form 

E ** "c [v f xB] (l-N f ) ' (14) 

v n n J 

As it can be seen from Eq.(14), in order to determine an electric 
field it is necessary to find the velocity of fluxoid motion v f 

for which purpose let us use the equation of motion suggested in 
Ref. [10] (when N n = 0, VT=0) : 

l-M ~ ^ m [ V f xk ] " ^ V s " *>» V f = °* ( 15 > 

where 77 ^ = n s e$ 0 /c, n g is a concentration of superconducting 
carriers; v g is their (average) current velocity; and 77 j_ and r? „ 
are viscosity coefficients. The existence of normal current i 

J n 

(for N n / 0 ) should have brought about the appearance of a 

v -dependent function A(v ) in Eq.(17). Since at N -» 1 all the 
n n n 

terms occurring in Eq.(17) go to a zero, the equation of fluxoid 

motion in this case come to relation A(v n )=0 coinciding with the 

equation of normal quasiparticles motion, i.e., to the Ohm's law 

(12) . We may assume that at N <1 variables v and v in the mo- 

n s n 

tion equation are also "divided", in other words, the equation is 
"split" into two expressions: A(v n )=0 (coincides with Eq.(ll) ) 

and Eq . (15) . 

Now let us discuss the way of modification of the fluxoid 
motion equation (15) with the presence of the temperature 
gradient. Firstly, "thermodynamic" force F , = -S,VT must 

T 1 0 

manifest itself in it. Here S^=S^(T) is a temperature dependent 

entropy of a fluxoid section of a unit length. This force is an 

analogy [2] for Lorenz force 77 [v xk] and can be derive from the 

'ms 

latter by substituting v g -» v*= (S^/r/ m ) [ VTxk] . Simultaneously force 

F t2 -analogous to "viscosity" force ~V±' v s from Eq.(15) should be 

introduced; it is derived from the latter by means of the same 
substitution v g -* v* : F t2 = -S^ (v±/V m ) [ VTxk] . 

Another reason for appearance of "thermal" forces in the 
fluxoid motion equation is the origination of counterflow 
supercurrent in-between their normal cores [11]. If there is a 
temperature gradient inside the fluxoid core (that is in the 
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region where N =1) there appears additional normal current with 
the density, according to Eq.(ll), equal to j T = - ^ N1 VT ~ 

-/? n2 [ VTxk] . This current is balanced with externally flowing 
counterflow current with an average density 3 b = j bs + j bn = -j T . 
which is a sum of supercurrent j fas and normal current j fan . The 
superconducting component of this current equals j fas = N g j b = 
-N j T . Similarly to the transport supercurrent j s = en s v s it is a 
source of forces effecting fluxoids [12]. Hence two more forces 

F T3 = -<V Ne ><0„,[VTxk]-0 N2 VT), (16) 

F„ = ( >lj./Ne)(^ |)1 vt+^ 2 [ VTxk ) 

should be introduced into fluxoid motion equation (15) related to 
the case when VT=0. These forces are derived from Lorenz force 
Tj m (v s xk] and viscous force ~V± v s by means of substitution v g -» 

-v**, where v**= j bs /en s = - ( 1 /eN) . 

Thus, in the presence of VT the fluxoid motion equation 
acquires the form: 


"I s * ■ (fc) < ' , j-^i + ’ , «^2 ) ] 7T ■ ".h xk ] ■ ,i ,v f = °' <17) 

Using expression j =en g v s = j- j (where the normal current follows 

Eq.(ll) ) and substituting relationship v g = c[ ExB ] [ ( l~N n f n ) /f ] /B 2 

ensuing from Eq.(14) into Eq.(17) we obtain the equation for the 
electric field. The solution of that equation in respect to E 
gives the relationship determining the transport coefficients of 
superconductors in the mixed state: a Seebeck coefficient 

(thermal power) 


j 

The negative sign appeared in the v g -> -v** substitution is a 

result of the reverse (in respect to the normal current ) 
direction of supercurrent j . The nondiscrepancy of the fluxoid 

motion equation (17) thus obtained is confirmed by the 
fulfillment of the Onsager relationship for calculated kinetic 
coefficients (see below). 
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a=N [a (o cr N +p a N )+Q B(p a* - p a N )"1 + (cS./S> )a (18) 

nL N ' r xx xx r xy xy' N Vf xx xy r xy xx'J v <f> 0 ,r xy v ' 


and Nernst coefficient 

Q=N k(P * +P )-(«»/B) 
^ n [ N ' r xx xx r xy xy' N 7 


N N 

p o - p o 
r xx xy r xy xx 


■]- 


( cS , /B$. ) p 
v <f> 0 /r xx 


(19) 


It is easy to see that for normal metal (N n =l, =0) Eqs . (18) 

and (19) give a=a N and Q=Q N » and in the superconducting state (no 

fluxoid motion, P<P C and p xx = p xy = 0) , as it might be expected, 

a=Q-0. Thus, Eqs. (18) and (19) are true for a whole transition 
region of the resistive mixed state of a superconductor. 


3 . Heat conductivity . 

Peltier, Ettingshausen and Righi-Leduc Effects 


Heat flow q g , related to the entropy transfer with moving 

fluxoids, equals q s = S^Tn ef v f , where n gf is effective density of 

moving fluxoids which is different from their real density 
n f =B/$ 0 - Nevertheless, rewriting Eq.(14) in the form E=-(l/c)[v 

x B](n /n f ) we come to the relation 

q s =T(cS^/S 0 )[Exk), 

which is an analogy to Eq.(6) . As to the heat flow due to the 
normal quasiparticle motion, it is defined by the relationship 

%.*"» VV‘« B < te in 1 +N n (-«„VT + L N B(kxVT]) , (20) 

which accounts for « N and L N being proportional to current 
carrier concentration N while IT,, and the product E k are 

n N N N 

independent on it. 

Thus, the electron component of the heat flow in a 
superconductor equals 

cS 

q=q s +q n = T(-^)[Exk]+n N j n +E N /c N B[kxjJ+N n (-/c N VT+L N B[kxVT]) . (21) 

To deduce corresponding kinetic coefficients is possible by means 
of expressing q in terms of total transport current j. First let 
us find the connection between normal and total currents 
substituting Eq.(l) into Eq.(ll): 


j =N [j(/> o N +p a N )+[kxj](/j <7 N - p a H )1- 

J n n [_ r xx xx r xy xy r xy xx r xx xy J 
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( 22 ) 


-(£ N 1 -W VT - (/WWi™*] 


where ^ = ( a xx a+a xy QB) » 0 2 = {a xy a ' a xx QB) • Next » substituting 

Eqs.(l) and (22) into (21) we come to expression (2) where 
Peltier coefficient 


IT = N [il (p a N +p a N ) +E„k..B ( p o" - p <j N )1 + T(cS ,/$> )p , , 
n [_ N v ^xx xx ^xy xy' N N VA xx xy r xy xx'J <p 0 r xy r 

Ettingshausen coefficient 

(E/cB) =N [e k B (p a N +p a N )-JI(p cr N -p a N )]-T(cS ,/<!>) p , 
v ' n |_ N N vr xx xx r xy xy ' N vf xx xy r xy xx'J <p 0 r xx ' 

heat conductivity coefficient 
and Righi-Leduc coefficient 

LB = N n L N B+n H^N2- N n^)- E N% B ^Nl- N n^)- T ( CS /°0^- 


(23) 


(24) 


(25) 


(26) 


Proceeding from the Onsager equation for a normal metal 
(IT =a,T, E„,/c =Q ,T) we make sure that Onsager equations (3) are 

also valid for above given Seebeck and Peltier coefficients 
(compare Eqs.(18) and (23)) as well as for Nernst and 

Ettingshausen coefficients (compare Eqs.(19) and (24)) for a 
superconductor in the mixed state. This is an evidence of 
inherent nondiscrepancy of the model suggested and, in 

particular, of a correct form of the fluxoid motion equation 
(17). The latter might be obtained without defining concretely 
expressions for "thermal" forces F T1 -F T4 • It would be sufficient 

first to obtain relations (18) and (19) for a and Q, proceeding 
from Onsager relations and Eqs . (23) and (24) for n and E 
coefficients, and then to "restore" the fluxoid motion equation 
(17) . 

The Hall conductivity (resistance) is known to be 
sufficiently lower than the longitudinal conductivity 

(resistance) at any temperature [3,4]. It gives a possibility to 
simplify noticeably the expression obtained for kinetic 

coefficients : 

a=IT/T * N n a H (^ x /^)+(c S^/<V P„ y * < 27 > 


Q=(E/kT) s N n Q„(p xx /^)-(cS 0 /Bt o )^ x * - (cS^/BJ. 0 )p xx , (26) 
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( 30 ) 



The first two of them were derived in earlier works in a not 
quite correct way [ 1 , 2, 11 , 13]. 

Let us emphasize that the expressions derived for 
coefficients a, Q, IT, E, k and L are universal, i.e., independent 
on the choice of coefficients f and f . The latter determine 

s n 


only the connection between total conductivity of a 
superconductor in the mixed state and "partial conductivities", 
related to normal quasiparticles and fluxoids [14]. 

From the expressions obtained for kinetic coefficients a, Q, 
II, E, k, and L of a superconductor in the mixed state it is seen 
that each of them consists of two terms, one of which is 
proportional to concentration N n of normal quasiparticles and the 


other is to the entropy of fluxoids. These two parts may be 

conditionally called "normal" and "fluxoid" contributions^. Their 
role in different kinetic coefficients is not identical. 
Longitudinal effects (described with coefficients a, II, and /c) 
are defined basically with the "normal" contribution, while 
transverse effects (coefficients Q, E, and L) are governed by the 
"fluxoid" one. It is due to the following: 

1) Seebeck and Peltier coefficients a, and II : a relatively 

minute value of the fluxoid entropy (cS .p../a v ,$> ~ R « p /p ) 

(p r N N 0 r xx r xy 

and a absolutely small Hall angle (P xx /P x y y> 1) t the vortices are 


almost immovable along the temperature gradient and transfer 
small (in comparison with normal electrons) energy; 

2) Nernst and Ettingshausen coefficient: a relatively large 
value of the fluxoid entropy ( cS^p^ /a N <X > 0 ~ R » 1 ) and a 

relatively small Nernst coefficient Q ki in the normal state 

J N 


(Q N B/a N ~ 7 b « 1 ) ; fluxoid transfer high (as compared to normal 

electrons) energy in the direction perpendicular to the 
temperature gradient; 

3) Heat conductance and Righi-Leduc coefficients /c and L: all 
the above-said conditions and a relatively high electron heat 
. 2 

conductance in the normal state (k„o „ /T a, » 1) : the vortices are 

N r N c N ' T 


4 

The conventionality of this division can be seen if only from 
the fact that the "normal" contribution into kinetic coefficients 
is related to the longitudinal resistivity P xx t which, at least 

to a certain extent, is determined by the fluxoid motion. 




almost immovable along the temperature gradient ant transfer 
small (in comparison with normal electrons) energy. 

A good agreement with experimental data for HTSC of different 
compositions has been achieved for a different selection of 

interpolation factors of the form f s =N^, f n =N^ [14]. The values p 

= v =1 for YBa_ Cu,0_ c and p - v =3/2 for Bi-2223 were shown to 

provide the best fit.. A possible reason for these different 
values is the different degree of anisotropy of the two compounds 
in question. This problem, however, is to be further analyzed. 
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